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Resume. We introduce a restricted four body problem in a 2+2 configuration extending the 
classical Sitnikov problem to the Double Sitnikov problem. The secondary bodies are moving on 
the same perpendicular line to the plane where the primaries evolve, so almost every solution is 
a collision orbit. We extend the solutions beyond collisions with a symplectic regularization and 
study the set of energy surfaces that contain periodic orbits. 



1. Introduction 

We extend the classical Sitnikov problem [4] that is a special case of the restricted three body 
problem to a restricted four body problem in a 2+2 configuration. This configuration consists in 
two massive bodies evolving on keplerian orbits around their center of masses and two infinitesimal 
bodies evolving on the perpendicular line that cross the center of masses. The Double Sitnikov 
Problem consists in determining the evolution of the infinitesimal bodies under the Newtonian 
attraction of the massive bodies. Since the evolution of the secondaries is coUinear we are interested 
in solutions with elastic bouncing at collisions. In this note, we consider the circular restricted case 
of the double Sitnikov problem that is the integrable case. We study the periodic orbits on resonant 
tori. 

2. The Circular Double Sitnikov Problem 

We consider that primaries are evolving in circular keplerian orbits and each one has mass 
nil = "^2 = 5- The secondary masses are related hy v = c fi with c £]0, 1]. We choose the 
normalized masses a = and f3 = 1 — a. With these conditions the total energy of the system 
depends on a and fi as parameters as follows 

where 

(2) M = 



a 

(3 



The Hamiltonian vector field is obtained by ixn^ — dH. The configuration space will he V = 
{93 > 94} the flow 1^9(2;, t) — (pt{x) is not complete by the singularity due to collision. To avoid 
the singularity and to extend analytically the equations to the hyperplanc = (74 we perform 
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a symplectic regularization with the mapping p : T*M? — > T*M^ defined through the generating 
function of second type 

12 ^ 



W{Q, p) = P3 [qa + /3^) + Pi [Qi - 



Then the mapping p : (Q, P) ^ (q, p) is such that p*{J2i dpi/\dqi) = J2i dPihdQi and, therefore 



p e S'p(T*R^). Also, we consider the time rescahng 4^ ~ otPQl- The regularized Hamiltonian 



function is i = {H — h) o p this Hamiltonian function depends on a, h and p as parameters 

and is valid only in the energy level L = for each h fixed. If z = (Q3, Q^, P3, P4) we have explicitly 

1 
2 



L = ^ (a/3 P2Q2 ^ p2) _ 2a/3V 



2a 213 

+ ■ 



y^(2Q4 + /3Qi)' + l ^{2Qi^aQlf + 1 



Obtaining lim^^o Lh{z] a, p) = Lh{z] a, 0) and reversing the process, since a(3Q'l is not identically 
zero, we recover the Hamiltonian function without the term p ^Jt^^ ■ Finally, to extend the solutions 
in a continuous way beyond collisions (with elastic bouncing) is necessary that the linear momentum 
be conserved and it is possible if and only if c = 1 . We have the following 

Proposition 2.1. In the circular double Sitnikov problem if p = v then the flow ft{x) of the 
limiting case p can be extended to a complete flow in a natural way. 

Then the Hamiltonian system of the double Sitnikov problem is the triplet (T*M^, oj, H) where 

4 111 
Lu dpi A doi, and H = -IpP , , P=(P3,P4). 

3. Periodic Orbits 

Theorem 3.1. The action-angle coordinates for the circular double Sitnikov problem takes the 
form 

f2 1 
(3) J{K) = ^{2E{h)~K{hC)-^{2klM)) B,(t-h,) = —t(v,k)^B^^, 

where fi^ — ^-jj^j^^(2P(fci) — A'(fci)+H(2fc,^, fc^)) is the return time of the secondaries, ki = 
and doijor i = 3,4, are constants determined by the initial conditions. 

Theorem 3.2. The solutions for the circular double Sitnikov problem can be written as 

fe,P3,<Z4,P4)(t) = ( l_2fc3^,2(,3) ^2^/2fc3 ^(-3), Y3^fcT^,2^/2fc4 c(.4)j , 
where vi are functions of t obtained inverting the function 

^ = / 4(1 -2^^71(1/02)2''''" 

and s{vi) = sn{i'i{t),ki), c{i>i) = cn{i'i{t),ki), d{vi) = dn{vi{t), ki) are the sine, cosine, and delta 
amplitude Jacobi elliptic functions, and ki = '•^2+^' Jqj- ^ = 3^4. 
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It is possible to integrate the expression (4) with cUiptic functions and clhptic integrals to obtain 



(5) t = ^ 



1 - 2fc2) 



„, ■, „, , 9n , ^ sn{v)cn(v)dn(v) 

2E{iy)-iy + n{iy,2k^)-Ak^- ^ ' ^ ' ^ ' 



where C is an arbitrary constant of integration. In [2] the reader will find a nice and complete 
study of this function. 

Definition 3.3. We say that ip{t) is a periodic solution of period r with t > if tp{t + t) = ip{t) 
for all t R and there does not exist f G (0, r) such that ip(t + f) = ip{t), i.e., r is the minimum 
period. 

We will use the notation (p, q,n) = 1 to mean that the greatest common divisor is gdc{p, q,n) = I, 
in other words, that the three numbers have not common factors at the same time. 

Proposition 3.4. For every periodic solution of the double Sitnikov problem there exist 3-plets 
{p,q,n) e Z'^ such that {p,q,n) = 1, and p > and p > holds. The periods of these solutions 
are related to the partial energies by t 2pi: ~ qT(hi) = nT(/i2). 

Definition 3.5. We say that an energy surface S/,., = H^^(h.^,) accepts a periodic solution if there 
exists {p,q,n) e N"^ with the properties : PI) {p,q,n) = 1, and P2) p > P > such that 

= T^^ (^g^) +2^^^ i.^ir) denote the set of the constant energy surfaces that accept periodic 
orbits as 

371 = {Ysh = H-^{h)\h = K) . 

Theorem 3.6. In the circular double Sitnikov problem there exists a countable number of energy 
surfaces S G 9Jl that contains resonant tori foliated by periodic orbits. Moreover, the set of values 
/i* G H{T*R^) C M such that S^, G is dense in (-4,0) and have zero measure in R. 

The proof of the theorem is an immediate consequence of the two following lemmas 

Lemma 3.7. For each n G N the circular double Sitnikov problem has periodic solutions of period 
2n-K 

We will just exhibit at least one periodic solution of period t = 2iV7r. This is immediate from 
the fact that there exists such periodic solutions in the circular (classical) Sitnikov problem. 

Proof.- For any G N we can choose the combination p ^ N and q = n = 1 that produce 
{p, q) ~ 1 and {p, n) = 1 with p > and p > and the proposition 2.8 in [2] assures that there 

exists hi,h2 G (-2, 0) such that T{hi) = ^ and T{h2) = ^ then the hypersurfacc 4-/12) 
contains a torus foliated by a family of periodic orbits with period r = 27tN = T{hi) ~ T{h2) k 

Lemma 3.8. For each iV G N fixed, the circular double Sitnikov problem have a finite number of 
tori foliated by periodic orbits with period t — 2N'k. The number 

m^{N)+ '^(9); with vip)^p\{{^-\^ 

(j<2V2Ar "|p 

is an upper bound (although is not an optimum bound). ip(p) is the totient function or Euler's phi 
function. 

Proof.- For each p G N fixed there exist 3-plets {p,q,n) G N"^, where properties PI and P2 of 
Definition 3.5 holds. Therefore, we search for the number Cp of 3-plets {p, q,n) = 1 coprimes. It 
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is easy to see that for every q < 2\/2p and {p,n) ~ 1, the 3-plet {p,q,ri) docs not have common 
divisors. This triplets are exactly {2V^p) ■ {2V2ip{p)) 8pip{p). 

Additionally, we must add all the couples {q, n) coprime such that {p, q) and (p, n) are not 
coprime. This means that for each integer q < 2\/2p with (p, q) ^ 1 we must add the number of 
coprimes ^p{q). Then we have 

(6) Cp<8pip{p)+ J2 "P^l) 

Finally we must eliminate the elements that are in both sets, however the number (6) is an upper 
bound of the triplets {p,q,n) £ where properties PI and P2 holds. 

The 3-plet {p,q,n) e induce a point x = {2n^,2TT^J e {T{h:i),T{hi)) such that the lagran- 
gian torus T= (/i-ioT"i)(a:) is foliated by periodic orbits of period 2A^7r, therefore it is a resonant 
torus Tijes C r*M2. A 

Proof.- (of Theorem 3.6) The first part of theorem is a consequence of the fact that the 
countable union of finite sets is a countable set. Using Lemma 2 and Lemma 3 we have that the 
number of resonant tori are countable, and since each torus belongs to exactly one energy surface, 
the set OJl is countable too. 

We define the map T : g* ^ by [h^^hi) ^ (^^^2^' ^^2^^) ■ '^''^^^ rational point y G 
Img{T) with y = (^, (r, s) = 1 and {u,v) = 1, we construct the point (y^, y-, y) S where 
g = gcd(ru, su,rv). Since this point fills properties PI and P2 of Definition 3.5, there exists a 
resonant torus foliated by periodic orbits with period r = 2-y7r = ^T{h^) = ^T{h4) The set of 
rational values of T defined by RP :~ Img{T)r\Q^ is a dense subset of zero measure in Img(T). The 
mapping T is continuous and then T^^(RP) C g* is a dense subset in the image of the momentum 
map p,. Now we construct the function 7i : g* ^ K such that sends x = (/13, 1-^/13 + /14. It 
is immediate that Ti{T~^{RP)) C (—4,0) is a dense subset by continuity, and have zero measure 
since RP is a countable set. ▲ 
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